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Abstract 

In this paper, we provide a theoretical description, and calculate, the nonlinear frequency shift, 
group velocity and collionless damping rate, v, of a driven electron plasma wave (EPW). All these 
quantities, whose physical content will be discussed, are identified as terms of an envelope equation 
allowing one to predict how efficiently an EPW may be externally driven. This envelope equation 
is derived directly from Gauss law and from the investigation of the nonlinear electron motion, 
provided that the time and space rates of variation of the EPW amplitude, E p , are small compared 
to the plasma frequency or the inverse of the Debye length, v arises within the EPW envelope 
equation as more complicated an operator than a plain damping rate, and may only be viewed as 
such because {uE p )/E p remains nearly constant before abruptly dropping to zero. We provide a 
practical analytic formula for v and show, without resorting to complex contour deformation, that 
in the limit E p — > 0, v is nothing but the Landau damping rate. We then term u the "nonlinear 
Landau damping rate" of the driven plasma wave. As for the nonlinear frequency shift of the EPW, 
it is also derived theoretically and found to assume values significantly different from previously 
published ones, assuming that the wave is freely propagating. Moreover, we find no limitation in 
kXu, k being the plasma wavenumber and Ad the Debye length, for a solution to the dispertion 
relation to exist, and want to stress here the importance of specifying how an EPW is generated to 
discuss its properties. Our theoretical predictions are in excellent agreement with results inferred 
from Vlasov simulations of stimulated Raman scattering (SRS), and an application of our theory 
to the study of SRS is presented. 
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I. INTRODUCTION 



Landau damping is a linear, collisionless process, resulting from the global acceleration 
of electrons by an electrostatic wave. Indeed, in the linear regime, an electron plasma wave 
(EPW) with phase velocity globally accelerates the electrons of initial velocity v < v^, 
and decelerates the other ones. When this leads to an overall acceleration of the electrons 
by the wave, as for example in an initially Maxwellian plasma then, because of energy 
(or momentum) conservation, the plasma wave damps while, in opposite regime when the 
electrons are globally decelerated, the wave grows unstable. The damping, or growth rate, 
h>L, of the EPW in the linear regime was first derived by Landau in his famous 1946 paper 
Ref . l|] . While addressing the growth of the EPW was rather straightforward, Landau had 
to use complex contour deformation and analytic continuation to derive the damping rate, 
a technique which initially shed some doubts into plasma physicists minds as regards the 
validity of Landau's calculation, all the more as Landau never clearly discussed the physics 
of the damping. Landau damping, or growth, is predominantly due to the nearly resonant 
electrons, those whose initial velocity Vq is such that |t> — < Vh/k, where k is the plasma 
wave number (while the exactly resonant ones, such that Vq = v^, do not contribute to it). 
Then, as is well known, if -C u) pe , where u pe is the plasma frequency, vl is approximately 
proportional to the derivative, /o(f</,), of the electron distribution function in the limit of a 
vanishing field amplitude. 

A nonlinear counterpart of vl was first calculated by O'Neil in Ref. pj, who considered an 
electron plasma wave of constant and uniform amplitude, E , which grew infinitely quickly 



in an initially Maxwellian plasma. When uj b ^> v Ll where ujb = \/ekEo/m, — e being 
the electron charge and m its mass, most of the nearly resonant electrons are trapped and 
oscillate in the wave trough. Within one oscillation period, a trapped electron neither gains 
nor loses energy in the wave frame, so that the mechanism which gave rise to Landau 
damping vanishes, and so does the damping rate after a few oscillations at a frequency close 
to u>b, as shown by O'Neil. 

A countless number of papers, addressing both the linear and nonlinear regimes, have 

n the linear regime, the physics 
4j and references therein), and 
new derivations of Landau damping which did not resort to complex contour deformation, 



been written since these two seminal works were publishec. 

n 

of Landau damping was extensively discussed (see Ref. [3|, 
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or which extended Landau's result to non smooth initial distribution functions fas is the 



case for a real plasma made of discrete particles) were found (see Ref. [3J]). Moreover, 
very recently, Belmont et al. showed in Ref. J5] the very unexpected result that an EPW 
could damp at a rate different, and lower, than that derived by Landau, provided that this 
wave was excited from noise in such a way that the electron distribution function had a 
complex pole in velocity space. This shows the importance of specifying the way an EPW 
has been created in order to correctly discuss its physics properties and, in particular, to 
correctly calculate its complex frequency. In this paper, we provide a derivation of the 
Landau damping rate which does not resort on complex contour deformation and which, we 
believe, is quite simple. This moreover allows us to discuss the ability to excite a plasma 
wave in such a way that it decays at a non-Landau rate. 

In the nonlinear regime, several papers recently discussed the very work of O'Neil, even- 
tually leading to its experimental check (see Ref. jj| and references therein). Although the 
situation considered by O'Neil is physical and could be reproduced experimentally, it is not 
the most general one since a plasma wave amplitude usually depends on both space and 
time, even when this wave induces nonlinear electron motion. Generalizing O'Neil's results 
has been a long standing problem in plasma physics, which we address in this paper. In par- 
ticular, we provide an analytic expression, supported by numerical results, for the nonlinear 
collisionless damping rate, u, of a plasma wave whose amplitude may vary in space and time, 
in the limit of non relativistic electron motion and slow amplitude variations. We moreover 
restrict to a driven plasma wave for the following reasons. First, only if an EPW is externally 
driven may it grow in an initially Maxwellian plasma and may global electron acceleration, 
at the origin of Landau damping, occur. Second, for a driven wave, the initial conditions 
can be defined unambiguously and, in particular, one may assume that the plasma wave 
amplitude is initially at a noise level. This allows one to discuss the generality of previous 
results, regarding the nonlinear dispersion relation of an EPW, derived by assuming that the 
was was freely propagating. Third, our work directly applies to stimulated Raman scattering 
(SRS), which is studied as a tool for amplification of electromagnetic radiation, but which 
may also be detrimental for an inertial confinement device such as the Laser MegaJoule 
because it may induce the reflection of a substantial part of the incident laser energy. Now, 
recent numerical 8, 9] and experimental [h]] papers on SRS reported reflectivities far above 
what could be inferred from linear theory. This so-called "kinetic inflation" was attributed 
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to the nonlinear reduction of the Landau damping rate, although no theory, nor analytic 
formula, was available to support this assumption. The present paper addresses this issue 
and discusses in detail the derivation and physics of the very recent results, published in 



Ref. 



Ill 



There are several caveats in trying to define, and calculate, the nonlinear collisionless 
damping rate, v, of a driven wave. For example, one cannot that easily use energy conserva- 
tion as O'Neil did, nor even momentum conservation, to derive z/, because the electrons are 
accelerated by both, the drive and the plasma wave. It is usually argued that the plasma 
wave amplitude, E p , is much larger than that, Ed, of the laser drive, and this argument has 
been used by Yampolsky and Fisch in Ref. [12j to derive a set of equations from which v 
could be derived numerically, in case of a purely time growing wave. The relative values 



of E p and Ed has actually been investigated in detail in Ref. [13| where it has been shown 
that only in the nonlinear regime when v rs 0, or in the linear regime when the Landau 
damping rate is small enough, is E p S> Ed- Moreover, even in these regimes, only the space 
integrated energy, or momentum, is conserved, and these global quantities are not easily 
related to v which is defined locally. Since v is not easily calculated using conservation laws, 
in this paper, we will derive it from Gauss law, which is unambiguous. Using the electron 
susceptibility, Xi introduced in Ref. [jjj], and whose definition will be recalled in Section HT| 
we find that, provided that when Re(x) ~ —1 and |Im(x)| <C 1 (which are easily achieved 
conditions), E p is related to Ed and to the dephasing 5ip between the plasma wave and the 
external drive by the equation, 

lm(x)E p - k~ l d x E p = E d cos(8(p). (1) 

Eq. ([I]) tells us how efficiently an electron plasma wave may be driven, which is an important 
issue since our work was primarily motivated by the estimating of Raman reflectivity in 
fusion devices. To this respect, the nonlinear derivation of lm(x), which will be discussed 
in detail throughout this paper, is essential since it is clear, from Eq. (pQ), that a nonlinear 
decrease of Im(x) would enhance the driving of the EPW and, hence, SRS. Now, it is also 
clear that, while it is driven, an EPW accelerates the plasma electrons exactly the same way 
as if it were freely propagating, which hampers its growth. The effectiveness of the EPW 
drive therefore significantly depends on the rate of energy, or momentum, transfer from the 
wave to the electrons, a process akin to that giving rise to the Landau damping of a freely 
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propagating wave. We would like to make this more transparent by writing Eq. (pQ) in terms 
of an envelope equation of the form, 



d t E p + v g d x E p + uE p = E d cos(^)/^Xenv ( 2 ) 

Then, v g would be called the group velocity of the plasma wave, and v its nonlinear Landau 
damping rate. In this paper, we indeed show how to derive Eq. (jSJ) from Eq. ([I]) and we 
actually provide an analytic formula for u, that matches the Landau damping rate, Ul, in 
the limit of vanishing field amplitudes. We moreover show that u, which depends on both 
the wave amplitude and its space and time variations, may be viewed as a plain damping 
rate because it assumes nearly constant values before abruptly dropping to zero. Then, 
not only is Eq. ([2]) physically more transparent than Eq. (pQ) but it is also easier to solve 
numerically to get, for example, quantitative estimates for Raman reflectivity. It is however 
important to note that the physical meanings of v and v g are not as obvious as for a freely 
propagating wave. Indeed, usually, the maximum of a driven plasma wave packet does not 
travel at v g . Moreover, the amplitude of the driven EPW does not decrease at rate u, but 
grows most of the time. Moreover, although Gauss law is unambiguous, there is actually 
not a unique way to write Eq. ([1]) in the form Eq. (j2j). However, because the transition to 
the regime where v « is quite abrupt, there is actually very little freedom in the choice 
for u, v g and d^xlnv m Eq. (EJ), which vindicates the use of that equation and the values we 
derive for its coefficients. 

The present paper, which is mainly devoted to the derivation of Im(x) and of the envelope 
equation Eq. (T5]), is organized as follows. For pedagogical reasons, we will first present in 
Section [Til the derivation of Im(x) in case of a purely time growing wave, and will explain 
how z/, v g and d^xlrw can be deduced from Im(%). In Section IHfl we will explain how 
these results can be generalized to a wave whose amplitude either grows or decays in time. 
Section [TV] addresses the issue of a time and space varying wave amplitude, and shows 
comparisons between our theoretical predictions and results inferred from one dimensional 
(1-D) simulations of SRS. In this Section will also be discussed how (3-D) effects may affect 
the range of the validity of the linear regime in terms on the EPW amplitude. In Section IVl 



we briefly recall results from Ref. [13J on the nonlinear frequency shift of a driven plasma 
wave, from which the dephasing 5(p stems and, in Section IVIl we show one example of 
the application of our theory to stimulated Raman scattering. Section IVHI concludes and 
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summarizes this paper. 



II. ENVELOPE EQUATION AND NONLINEAR LANDAU DAMPING RATE 
FOR A TIME GROWING DRIVEN PLASMA WAVE 

In this Section, we derive the envelope equation for an EPW whose amplitude only 
depends on time, and grows with time. This will allow us to introduce in a simple way 
most of the concepts useful in the general situation of a time and space dependence of the 
wave amplitude. Most of this Section is devoted to the derivation of Im(x), performed by 
using two very different methods yielding values of Im(%) which do match over a finite range 
of wave amplitudes. For small amplitudes, we use a perturbative analysis which provides 
an expression for lm(x) that clearly shows how v decreases as more and more electrons 
are getting trapped in the wave trough. Then, when v « 0, one can approximate Im(x) 
by, Im(x) = FpduXlnvi where T p is the wave growth rate, T p = E~ l dE p / dt, and c^Xenv 
is calculated by making use of the adiabatic approximation. As is illustrated in Fig. [2, 
the "adiabatic" and perturbative estimates of Im(x) assume very close values over a finite 
range of wave amplitudes, which allows us to derive an expression for Im(x) valid whatever 
the wave amplitude by "connecting" the two previous estimates, as shown in Fig. HI This 
connecting is made through a Heaviside like function, leading to abrupt changes in the 
coefficients of the envelope equation, Eq. (J2]). In particular, v is found to assume nearly 
constant values before abruptly dropping to 0, and this drop is concomitant with a sudden 
rise in d^xlnv ( see Fig. E]). Indeed, as will be shown here, that part of Im(x) which, in the 
linear regime, provides u, renormalizes c^Xenv when v ~ 0. 

Let us now enter the details of the theory. We consider here a driven plasma wave, 
meaning that the total longitudinal field (along the direction of the wave propagation) is the 
sum of the EPW field, which is a genuine electrostatic field induced by charge separation, 
and of the driving field (the so-called ponderomotive field in case of laser drive). We assume 
that both the electrostatic, E e i(x,t), and the driving, Ed r ive(x,t), fields can be expressed in 
terms of a slowly varying envelope and an eikonal, that is, 



E d (x,t) = E p (t)sm[(p p (x,t)), 
E drive (x,t) = E d {t) cos[v9 p (x,t) + 5(p(x,t)] 



(3) 
(4) 
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with \E~ld x E Pid \ < \k\, k = d x (p p , \E~^ d d t E p4 \ < |a>|, to = -d t (p p , and \5cp\ < \(p p \. Then, 
the total longitudinal electric field, including the plasma wave and the drive, also writes in 
terms of a slowly varying envelope and an eikonal, 

E d + E dnve = EoWeW'ti + c.c, (5) 



where Eq and <p are given in terms of E p , E d , (p p and 5p in Ref. 141 ] . This total field induces 
a charge density which may therefore be written as, 

p(x,t) = p {t)e^ + c.c. (6) 

Throughout this paper we assume immobile ions, and define the electron susceptibility as, 

X = T^Eq (7) 

When the plasma wave is not driven and Eq is an electrostatic field, then Gauss law straight- 
forwardly yields the usual dispersion relation 1 + x = 0. In the general case, we use the total 
field amplitude Eq in the definition of x so that the expression of the electron susceptibility 
would be the same, in terms of the field amplitude and of the unperturbed distribution 
function, whether the wave is driven or not. In particular, it is easy to show that, in the 



linear limit, x is nothing but the usual linear electron susceptibility, as derived in Ref. 15]. 
Plugging Eq. ([7]) into Gauss law one easily finds, 

lm{ X )E p = E d cos{5if), (8) 

provided that Re(x) ~ —1 and |Im(x)| 1 (see Ref. jl^] for details). In order derive 
Im(x) and cast Eq. ([8]) in the form of the envelope equation, 

d t E p + uE p = E d cos(fy>)/d w Xmv> (9) 

we now need to express x m terms of the electron distribution function. From Eq. ([TJ) , it is 
clear that po is nothing but a Fourier component of p so that, 



-ne 



f(tp,v,t)e tip dvd(p 



IX J — OO 



2tt 

-ne(e^) (10) 
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where n is electron density, / is the electron distribution function normalized to unity, and 
(.) stands for a local, in space, statistical averaging. For the sake defmiteness, and without 
loss of generality, we henceforth assume that E Q is a pure imaginary number, so that Im(x) 
is proportional to (sin(<£>)). As a first step to calculate (sin(<£>)), we need to evaluate which 
electrons significantly contribute to it. This is done by investigating the electrons orbits in 
phase space, schematically displayed in Fig. [IJ If E$ were a constant, these orbits would 
be exactly symmetric with respect to the velocity axis, and (sin(yj)) would be 0. Since Eq 
slowly varies with time, the electrons orbits are slightly non symmetric, all the more as the 
growth rate of the total field, T = E^dEo/dt, is small compared to the time it takes for ip, 
or the polar angle in phase space, to change by 2tt. This time, henceforth termed the pseudo 
period of the orbit, is very close to 2tt/ujb for a trapped orbit far enough from the virtual 
separatrix (which is defined by freezing the wave amplitude). Hence, as shown in Fig. [IJ 
when wg > T deeply trapped orbits are nearly symmetric with respect to the f-axis, and 
electrons on such orbits contribute very little to (sin(y)), and therefore to Im(x). This lets 
us derive a specific criterion as regards the electrons which need to be accounted for when 
calculating Im(x). Since Eq varies slowly with time, we use the adiabatic approximation 
to find out which electrons are trapped in the wave trough. In terms of the dimensionless 
wave amplitude $ = eEo/kT e , where T e is the electron temperature, and of the electron 
initial velocity, vq, and wave phase velocity, v^, both normalized to the thermal velocity, 
Vth = \jT e /m, the condition for trapping derived from the adiabatic approximation is, 
|fo — Vtj\ < 4v^$/vr (see Ref. jl4|). Then, an electron orbit will be considered as "deeply 
trapped" if \vq — < 4v / $/vr(l — SV), with SV large enough for the electron orbit to be 
nearly symmetric. Since the symmetry of a trapped orbit is governed by T/ujb, we choose SV 
proportional to 7/v^ = T/lub where, in order to stick to dimensionless variables, we have 
defined 7 = T/kv t h- Therefore, we will henceforth assume that an orbit is deeply trapped, 
and that the electrons lying on it contribute very little to Im(x), if |i> — v<f\ < V, with 
Vi = max jo, (4y / $/7r) 1 — 37/2 j, where the value 3/2 has been found numerically (see 
Ref. [jjj]). It is noteworthy that, for a slowly growing wave, 7/v^ = T/uob ~ 2/ j usdt, 
so that Vi > when j UBdt > tt, that is after the first trapped electrons have completed 
about one half of their pseudo periodic orbit and the phase mixing process, introduced by 
O'Neil in Ref. to explain the nonlinear decrease of v, has started to be effective. 

Let us now explain how we actually calculate Im(x) from the matching of two different 



S 
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FIG. 1: Orbits of electrons acted upon by the longitudinal field 2Eq sin^), whose amplitude slowly 
varies with time. The dashed curve is the virtual separatrix. 
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FIG. 2: — (sm((p))/& as a function of \f$> calculated numerically (blue solid line), pertubatively 
(green dashed line), and adiabatically (red dashed-dotted line) when the normalized wave phase 
velocity is = 3 and the normalized growth rate is 7 = 0.01. 

estimates. For small wave amplitudes (and more precisely when C 7), we use a per- 
turbative analysis to derive Im(x), while when \f$> ^> 7 (or more precisely Vi ^> 7) we will 
show that Im(x) is nearly proportional to 7 and can be very accurately estimated by using 



9 




FIG. 3: — (sin(</j))/3> as a function of \/$ calculated numerically (blue solid line), pertubatively 
(green dashed line), and adiabatically (red dashed-dotted line), when the normalized wave phase 
velocity is v $ = 4 and, panel (a), when the normalized growth rate is 7 = 0.1, panel (b), when the 
normalized growth rate is 7 = 0.2. 

the adiabatic approximation. Let us start with the perturbative estimate of Im(x). There 
are several reasons to believe that a perturbative analysis will be useful in deriving Im(x). 
First, it has been proven in Ref. [16( that for small enough wave amplitudes linear theory, 
which stems from a first order perturbative analysis of the electron motion, is valid. Second, 
as shown before, one may neglect the contribution of the deeply trapped electrons, whose 
motions cannot be treated perturbatively, to estimate Im(x). Mathematically, this amounts 
to bounding from below the small denominators in the perturbative expression of Im(x). 
Actually, although rigourous estimates remain to be done, it appears from the results of Ref. 
fl^ that the "small parameter" of the perturbative expansion of Im(x) varies from v^/7 
when Vi <C 7, to y/^/Vi when Vi 3> 7, and hence remains bounded. However, a perturbative 
estimate of Im(x) eventually ceases to be accurate as the wave grows. Physically, this may 
be understood by the fact that, as v$/7 increases, the electrons have to lie on orbits closer 
to the separatrix to significantly contribute to Im(x), and the motion close to the separatrix 
is known to be non perturbative. Note, again, that v^/7 ~ / ujBclt/2, so that \/$/7 be 
large corresponds to the usual criterion for a highly nonlinear, and hence non perturbative, 
electron response. 

Let us now detail the perturbative expression of Im(x), which sheds a lot of light on the 
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(b) 11 th order 




2 4 6 8 1( 



FIG. 4: — (sin.(</?)}/$ as a function of v^/7 when the normalized wave phase velocity is v s = 3 and 
the normalized growth rate is 7 = 0.01, as calculated numerically (green solid line) and theoretically 
using Eq. (fl"6j) (blue solid line) with, panel (a), Im(x pe r) calculated using a 1 st order perturbation 
analysis, panel (b), Im(x pe r) derived from an 11 th order perturbation theory. 



nonlinear decrease of z/, and on how Im(x) may be estimated when v$/7 is large. Indeed, 
at first order in the perturbation analysis, and at 0-order in the time variations of 7 and v^, 
one finds the well known result (see Ref. [yj] for details), 

-2 f 7t> 



Im(x) 



{k\ D f J M>Vl ( 7 2 + v*y 



fo(v + v^dv, 



(11) 



where \d = v t h/u pe is the Debye length and fo is the normalized electron distribution 
function in the limit of a vanishing field amplitude. Now, Eq. ([8]) derived from Gauss law 
is the envelope equation Eq. ([9]) only if Im(x) may be written as, Im(%) rs 5Ii + T p 5I 2 , 
where T p = E^dEp/dt, and where 8I\ and SI2 only depend on the wave amplitude and not 
on its time variations (at least over finite ranges of amplitudes). As shown in Ref. 13), 
E p 3> Ed, so that E p « E , except if v is so large that the term vE v dominates in the left 
hand side of Eq. ([9]) and E p « E&jv. Hence, T p w T = E 1 dE /dt. We therefore only 
need to write Im(%) as Im(x) « + r<5/2, which clearly requires to isolate the divergence 
of the integrand in Eq. (fTTj) when VJ = and 7 — > 0. We do this by using the following 
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FIG. 5: Panel (a), variations of d^xlnv an d, panel (b), variations of u, normalized to their linear 
values, as a function of y/Q/j, calculated when = 3 and 7 = 0.01 by using a first order 
perturbative analysis to derive Im(x p0 r)- 



decomposition, Im(x) = l\ + J 2 , with, 

r _ -2/oK) 
j-i = 



'JV 



(k\ D ) 2 i H >y i (7 2 + ^) ; 
fo(v 4 



-dv 



(k\ D ) 



7 



7 2 + ^j 



/9 = 



-27 



:[/o(v + ity) - vf' Q {v^)]dv. 



(kXn) 2 y w> vs (7 2 + ^ 2 ) 2 ' 
Since 7 -C 1, one may approximate / 2 by replacing ( 7 2 + t> 2 ) by t> 2 to find, 



-2 



r7 



(k\ D ) j\ v \>Vi 
= T(dxl/du), 



dv 



(12) 
(13) 



(14) 
(15) 



where the integral in Eq. (fT41) has to be taken in the sense of Cauchy's principal part when 
^ = 0. 

When Vi = 0, Xi is just the adiabatic approximation of the linear value of Re(%) (and 
its value does not change much provided that Vi < 1), while Ii = —-niJiX^^f^v^). Hence, 
Im(x) is in the desired form, Im(%) = — -n{k,\£>)~ 2 f^v^) + r<9 w Xi\ so that Eq. (jSj) may indeed 
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be written as the envelope equation, Eq. (jHJ), with d^Xem = ^Xi> an d v = u l, the Landau 
damping rate in the limit vl "C oj pe . 

When Vi 3> 7, I\ is nearly proportional to 7 and therefore so is Im(x), which is actually 
obvious from Eq. ffTTT) . Then, Eq. straightforwardly writes as Eq. with v « 0; 
Landau damping has become negligible in the time evolution of the driven plasma wave. 
Physically, the decrease of v towards is due to the trapping of the nearly resonant electrons, 
which no longer contribute to v while oscillating in the wave trough, just like in the situation 
considered by O'Neil. 

Replacing ( r y 2 + v 2 ) by v 2 in Eq. ( fTTj) . which is valid when V\ ^> 7, one actually finds 
Im(x) = r^Xefr) 1 ) where (Xcfr) 1 ^ s the first order estimate of some effective real suscep- 
tibility, calculated adiabatically and by removing the contribution of the deeply trapped 
electrons. There is however no need to resort to perturbation analysis to calculate xls s i nce 
this can be done by directly using the adiabatic approximation, as shown in Ref. [ijj]. 
Then, from the previous discussion, we expect that when a/$ S> 7, Im(x) « Td^xls, which 
provides a non perturbative estimate of Im(x), which will henceforth be termed the "adia- 
batic estimate" of Im(x) [although this is not a proper terminology since a direct adiabatic 
calculation of Im(x) would just yield Im(%) = 0]. It is noteworthy that the I\ term Eq. 
( TT2|) which, in the linear limit provides z/, fully contributes to d^xls when v « 0. In the 
strong damping limit, when 3> T, d^xlm m ay then increase by more than one order of 
magnitude, as illustrated in Fig. 

Let us now compare the perturbative and adiabatic estimates of Im(x) to those derived 
from test particles simulations. Numerically, we calculate the dynamics of electrons acted 
upon by an exponentially growing wave and estimate (sin(y?)) = J2f=i w i ^(Vi)) where the 
sum runs over all the electrons used in the simulation, and Wi = fo(voi), where vo% is the 
initial velocity of the i th electron and fo is the normalized unperturbed distribution function. 
In our simulations, we chose fo(v) = (27r)~ 1//2 exp(— v 2 /2). Whatever the wave phase velocity 
and for small enough growth rates, we always found that the high (11 th ) order perturbative 
estimate of Im(x) was valid at least up to \/¥/7 ~ 10, while the adiabatic estimate was 
correct whenever \/$/7 > 3 (see Fig. [2]). Such comparisons moreover allowed us to conclude 
that an adiabatic estimate of Im(x) was only accurate if 7 < 0.1, as illustrated in Fig. [3j 

Using the perturbative, Im(x per ), and adiabatic estimates of Im(x) within their respective 
ranges of validity, which do overlap, we obtain the following expression for Im(x), valid 
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whatever the wave amplitude, 

Im( X ) = Im( Xper ) [l - Y (V$/ 7 )] + rd uX r eS Y (v^/ 7 ) , (16) 

where Y is a function rising from to 1 as \/$/7 increases. Since, as shown in Fig. [2j the 
convergence of Td^xls towards Im(x) is quite sharp, Y should rise very quickly from to 
1 as v$/ 7 increases from a little less than 3 to a little more than 3. This is the case if we 
choose Y(x) = tanh 5 [(e I//3 — l) 3 ]. Fig. H] shows comparisons between theoretical values of 
— (sin (<£>))/$ derived from Eq. ffTBT) . and numerical ones provided by test particles simula- 
tions. From this Figure, it is clear that using a high (11 th ) order perturbative expression for 
Im(x per ) yields very accurate values for — (sin((^))/$, and hence for Im(x), while calculating 
Im(x pcr ) at first order already yields very good results, with much more simple formulas! 
Therefore, for practical purposes such as the numerical study of SRS, we restrict to first 
order expressions. Then, from Eq. ( |T6l) and the expression found previously for Im(x per ), 
we conclude that Gauss equation, Eq. (jSJ), is the envelope equation, Eq. (jUJ), with, 

X r em = (l-Y)x x [ + Yx X : s , (17) 

V = Y X ARXenv « Y X h/d^l, (18) 

where I\ and Xi are defined by Eqs. ( |T2i) and ( |T5l) . If we were to replace 7 by 
(kvthE p )~ 1 dE p /dt in the expression (|T2l) for Ji, we would find that actually is much more 
complicated an operator than a plain damping rate. However, as shown in Fig. El provided 
that 7 remains nearly constant, v assumes nearly constant values before abruptly dropping 
to 0. v may then indeed be viewed as a damping rate, both physically and when numeri- 
cally solving the envelope equation, Eq. ([9]). We therefore successfully defined an effective 
nonlinear damping rate, u, yielding the time evolution of the driven plasma wave, which 
was our prime goal. We term v the "nonlinear Landau damping rate" of the driven plasma 
wave because it physically stems from the electron acceleration by the EPW, which is the 
very mechanism giving rise to the Landau damping of a freely propagating wave. Then, as 
expected, the linear value of v is nothing but the Landau damping rate. Note that we relate 
v to the growth of the driven plasma wave and not to any other quantity, such as the energy 
gain by the electrons from the wave. As shown in Fig. El the drop in v is concomitant with 
a rapid growth of d w x r e mi since the term in Im(%) which gives rise to v in the linear regime 
fully contributes to d^xlnv when v ~ 0, so that Im(x), and the efficiency of the driving of 
the EPW, varies smoothly. 
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III. GENERALIZATION TO AN ARBITRARY TIME DEPENDENCE OF THE 
WAVE AMPLITUDE 



In this Section, we generalize the results derived previously to a plasma wave whose 
amplitude may vary arbitrarily in time, provided that the growth rate, T = E^dEo/dt, is 
still such that |T| <C uj pe . We shall moreover show that the formula f|T8|) for z/, with I\ given 
by Eq. (|T2|) . is still useful provided that 7 is defined properly i.e., by Eq. (|33|) . 

We start by estimating (e~ lip ) through the means of a first order perturbation analysis, 
which proved in the preceding Section to be an important step in the derivation of Im(x). 
By using the Hamiltonian perturbation analysis detailed in Appendix \K\ one finds, at first 
order, ip(r) = ip + (v — v^)r + 5(p, where r = k\ D u pe t, velocities are still normalized to the 
thermal one, and, 

dw J 

where we have denoted w = v — v^. Then, 

= -[ fo(w + v^)^- [ Hu)e^ u -^dudw, (20) 

where V\ is a straightforward generalization, using the phase mixing argument, of the value 
found in the previous Section i.e., Vi = ^\f^ji: 1 — 3/ J l uj B (u)du , and where f is the 
electron distribution function in the limit $ —>■ 0. If $ has kept on increasing with time, fo 
is nothing but the unperturbed distribution function. If $ has reached a large enough value 
to induce nonlinear electron motion before decreasing to nearly 0, a perturbative analysis of 
the electron motion from t = is no longer valid. However, one may calculate the electron 
motion perturbatively from t = +00 by making use of the time-reversal invariance of the 
dynamics. Then, / is the distribution function in the limit t — > +00 which, as shown in 
Ref. [14J |. and as illustrated in Fig. [61 results from symmetric detrapping. In the interval 
\v — t^l > max(VJ), fo(v,t = +00) assumes the same values as the initial, unperturbed, 
distribution function, while in the interval \v — < max(V/), fo{v,t = +00) is nearly 
symmetric with respect to v^. Then, electrons whose initial velocity lies within the latest 
interval contribute very little to Im(x). This means that once deeply trapped, electrons no 
longer contribute significantly to Im(x), even after being detrapped. Eq. (!20l) may therefore 
be simplified by using for f the unperturbed distribution function and by replacing Vi by 
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max(Vz). Such a simplfication will be implicitly used throughout the remainder of this paper. 




FIG. 6: Results from Valsov simulations of stimulated Raman scattering showing, panel (a), the 
plasma wave amplitude (in its reference frame and in arbitrary units) as a function of time (normal- 
ized to the laser period), panel (b), the space averaged electron distribution function at the three 
different times indicated by the arrows in panel (a). Note that, as the EPW amplitude decreases, 
the space averaged distribution function becomes more symmetric with respect to v^, and that it 
is not only a function of the EPW amplitude. 



We now use the same kind of decomposition as in the previous Section to find a suitable 
expression of Im(x) i.e., we write, x — ~ i{k\D)~ 2 (e~ ltp ) / ® = Xa + Xb, with, 

*» = TTT^if S f *(«)(« ~t) [ iwe^dwdu, (21) 
(A;A D ) 2 $(r) J Q J H>Vl 

(JL J $ w e ™ (u_T)rfu ) dw - ( 22 ) 

Provided that $(r) 3> $(0), integrating Eq. (|22|) by parts with respect to time yields, at 
first order in the time variations of $, 

Im( X6 ) - -2(kX D )-^f [ f^ + ^)- W f>(v,) dw 

dr J\ w \ >Vl w 6 
= Tidxl/du) (24) 

where, again, the integral in Eq. ( J24l) has to be taken in the sense of Cauchy's principal 
part when V\ = 0. Hence, in the limit of a slowly varying wave amplitude, the expression of 
lm(xb) is exactly the same as that of the term J 2 found in the previous Section, Eq. (fL4")) . 
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When Vi — 0, since u — t), where 8 is the Dirac distribution, 

one easily finds Im(xa) = — (^Ad) 2 ^/^^). Hence when V\ — 0, which corresponds to the 
linear limit, Im(x) = Td^xl ~ 7r (^^i)) _2 /o(' y </')' so ^ na ^ Gauss equation is, 

d t E p + v L E v = E d cos(^)/^Xi, (25) 

where vt is the Landau damping rate, in the limit vl <C u pe . Since our linear calculation is 
valid whether the plasma wave is driven, or not, it unambiguously shows Landau damping, 
without resorting to complex contour deformation. This therefore allows us to conclude that 
non-Landau damping, as described by Blemont et al. in Ref. cannot be obtained by 
using a drive at the same frequency as the plasma wave to excite it above the noise level, 
and then let it freely propagate. 

In the nonlinear regime, and when Vf l is much smaller than the typical timescale of 
variation of $, t$, calculating the time integral in Eq. (|2"T1) by parts yields, 

lm(xa) = -(k\ D r 2 ^ l f^)[Wf l d^/dT + 0(Vr 3 d 3 ^/dr% (26) 

Hence, when VJ ^> r^ 1 , Im(x ) is nearly proportional to T, and therefore so is Im(x), which 
implies v ~ 0. Again, as in the previous Section, we find that the decrease of v towards is 
due to the trapping of the nearly resonant electrons. Moreover, it is easy to show that in the 
limit Vi ^> 7, the I± term Eq. ( fl2l) of the previous Section, is close to — 4(k\]j)~ 2 f^v^^ /Vi, 
just as Im(x a ). We therefore conclude that, in the limit of large Vi, when !/«0, the results 
obtained in the previous Section for a growing wave, are valid whatever the time dependence 
of the wave amplitude. Hence, when Vi ^> r^ 1 which, for a slowly varying wave is typically 
the case when J^cosdu 1, we expect Im(x) ~ Td^xlsi where xls is the same as in the 
preceding Section. Then, generalizing the results of Section [Til we propose the following 
expression for Im(x), 



Im(x) = Im(xpcr) 



Y(2^u B dv) +Td ujX : s Y( K 2^uj B duy (27) 



where Y is the same function as for a growing wave, and when lm(x pcT ) is still the pertur- 
bative estimate of Im(x) which, at first order, is Im(x per ) = Im(x a + Xb) defined by Eqs. 
( I2TH221 . Eq. ( 1271) . when generalized to allow for the space variation of the wave amplitude, 
yields results in very good agreement with those inferred from Vlasov simulations of SRS, 
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as shown in Fig. [7J From the expression (l2Tj) of lm(x), and Gauss law Eq. ([8]), we derive 
the envelope equation (Q with, 

X r em = (l-Y)x x [ + Yx X : s , (28) 
v = Y x Im( Xa )/6Uonv- (29) 
Since, whenever Y is not close to 0, xLv ~ Xi> Eq. (I2"9"j) for i/ may be simplified in, 

!/ « y x Jm(xa)/d uX v (30) 

We now try to find a more simple expression for Im(x a ), leading to a practical analytic 
formula for v. In the limit of large Vi, we already showed that Im(x a ) was well approximated 
by Eq. ffl2l for I\. In the opposite limit when V\ <C r^ 1 , as shown in Appendix [Bj we find, 
Im(xa) = -(^A d )" 2 /o(^)[tt + Sxa], with, 

s ^-S)[fJl* {m '^- (3i) 

Similarly, when V; 7, a Taylor expansion of Eq. ( fl2l yields Ji = — (&;A£>)~ 2 /Q(f0)[7r + <5/i], 
with 5/i ~ — (4/3)(VJ/7) 3 . Since, for a slowly varying wave, and when $(r) ^> $(0), 
<5Xa ~ — [4VJ 3 /3$(r)] (J T $(n)c/u) 3 , we find that Eq. (TP2"|) still applies in the general case, 
and in the limit VJ <C r $ , provided that 7 be replaced by $(r)/ J r Hence, while for 

an exponential growing wave, for which Eq. ( fl2l is exact, 7 = $ _1 c?$ /dr = $(r)/ J T $(ii)cfw, 
we find that this equation still holds in the general case provided that, 7 = 3>(t)/ JJ" $(m)c?m 
when Vi <C r^ 1 , and 7 = $ _1 <i$/(ir when VJ 3> t7 . Therefore, we propose the following 
approximate expression for Im(x a ), 



Im(xc 



(A;A D ; 



2 



^/VA 2 7 ^ 
— 2ten UJ+T^ 2 



(32) 



_ $(r) - $(r - Tr/Vfl 
7 " ^ ' ^ 

where it is clear that 7 defined by Eq. (133]) has the required properties, 7 ~ $(t)/ /J" ${u)du 
when <C r^ 1 , and 7 ~ $ _1 d<I>/<ir when V; ^> t7 . Eqs. ( T32lT33|) have been used when 
comparing our theoretical estimate to numerical ones, and the good agreement between these 
two estimates, illustrated in Fig. [7J, shows the relevance of our approximation. Then, Eq. 
( 1301) . together with Eqs. ( 13211331) . provide a practical analytic formula for v. The accuracy 
for Im(x), and thus for i/, can even be improved by using, instead of Eq. (1321 . a result 



derived at higher order in the perturbative analysis (see Ref. ([14|) 
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IV. SPACE AND TIME VARIATION OF THE WAVE AMPLITUDE 



A. One dimensional (1-D) space variation and comparisons with 1-D simulations 
of Stimulated Raman Scattering 

The results of the previous Section are easily generalized to allow for one dimensional 
(1-D) space variations of the EPW amplitude. Indeed, using a Fourier expansion of the 
charge density, one finds, as shown in Ref. [yj], 

Im(x) = Im(xoo) - 4 d kX r ea y + Re(x)A], (34) 

where n = E^dxEo ps E~ l d x E p , and where Im(xoo) is given by Eq. (1271) except that all 
quantities must now be evaluated in the wave frame. More precisely, J^usdu in Eq. (I2T|) 
or in the definition of VJ now is, J*ljb[x — J^v^t^dt' \u]du, and the value for 7 to be used 
in Eq. ([32]) is, 



$(x, r) - $ 


x - 


- j T T _^ /Vi V^ U )du,T-Tl/Vi 


It-tt/Vi 


$ 


x ~ S u v 4>{ t ')d t 'i u 


du 



7(2,7-) = — p = -. (35) 

f^ /Vl $[x- Jl v^(t')dt>, du 

Plugging Eq. (1341) into Gauss equation ([T|), we find the following envelope equation, 

d t E p + v g d x E p + uE p = E d cos(<fyO/d w x£w> ^ 36 ) 

where, provided that [1 + Re(x)} w 0, v g = -d k x r cm / d^xl^ = uj/k - 2/[kd w xU\- lt is 
noteworthy that, since in the nonlinear regime xlnv R- e (x) ? v g dui/dk. Actually, since 
d^xlnv m ay reach values much larger than in the linear limit, the nonlinear value of v g may 
get quite close to the EPW phase velocity, as shown in Fig. [S] (d). Moreover, duj/dk may 
actually change sign from positive to negative, at rather small wave amplitudes (if k\n is 
large enough), which would entail a shock in the plasma wave profile if v g were indeed du/dk, 
while such a shock is not observed in kinetic simulations of Stimulated Raman Scattering. 
This is an indirect evidence that v g 7^ du/dk. 

We now compare our theoretical calculations of Im(x) against direct 1-D Vlasov simu- 
lations of SRS using the Eulerian code ELVIS |9j. In our numerical simulations, which are 
detailed in Refs. 0Jl3], the EPW results from the interaction of a pump laser, entering from 
vacuum on the left (x = 0), and a small-amplitude counterpropagating "seed" light wave 
injected on the right. Using a Hilbert transform of the fields, one can numerically calculate 
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FIG. 7: Time variations of Im(x) as calculated theoretically (green dashed line), and as calculated 
numerically (blue solid line) without, panel (a), or with, panel (b), using a Lorentzian factor in the 
v x B term of Vlasov equation. 

the ratio [E^ cos(<5<£>) + fc -1 t9 x ^7p]/ E p , which from Eq. ([T]) yields a first, numerical estimate, 
of Im(x). From Vlasov simulations one can also extract the values of all the quantities, 
such as J ujBdt and 7, which enter our theoretical formula for Im(%). Using these values we 
calculate a second, theoretical estimate, for Im(x). Both these estimates are compared in 
Fig. [7J plotting Im(x) as a function of uit, where ui is the laser frequency. The simulation 
results of Fig. |7] correspond to a plasma with electron temperature, T e = 5keV, and electron 
density n = 0.1n c , where n c is the critical density. The laser intensity is Ii = 4 x 10 15 W/cm 2 
while the seed intensity is I s = 10 _5 J/ and the seed wavelength is A s = 0.609/im. The 
results plotted on Fig. [3 (a) correspond to a simulation box of length L = 285 A/, where 
Az = 0.351/zm is the laser wavelength, and were measured at x — 77A/. In case of Fig. [3(b), 
the length of the simulation box is L = 350 A/, while the data were measured at x — 150A;. 
Moreover, in case of Fig. [3(b) , the v x B term in Vlasov equation was artificially multiplied 
by a Lorentzian factor, so as to mimic laser focusing which would occur in more than one 
space dimension. As can be seen in Fig. [7[ there is a very good agreement between the 
theoretical and numerical values of Im(x), especially as regards the decrease of Im(x) from 
its linear value in Fig. [3 (a), while the oscillations in Im(x) due to those of 7 are very well 
reproduced in Fig. [3 (b). 

The time variations of all the terms in the envelope equation ( |36|) are plotted in Fig. [8] for 
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FIG. 8: Panel (a), Im(x) as calculated numerically (blue solid line) and theoretically (green dashed 
line, panel (b) v normalized to the Landau damping rate, panel (c), c^Xenv normalized to its linear 
value and, panel (d), the EPW group velocity (blue solid line) and phase velocity (red dashed line) 
normalized to the thermal one. 

the same conditions as in Fig. [7] (a). Fig. [SJb) clearly shows that v remains nearly constant 
before abruptly dropping to 0, and that this is concomitant with a sudden rise in d^xlmn as 
for a purely time growing wave. This is very different from the oscillating result found by 
O'Neil because, in this paper, we consider slowly varying waves inducing a nearly adiabatic 
electron motion. As a consequence, electrons orbits are deformed as the wave grows so that 
electrons with the same initial velocity are all trapped nearly simultaneously, and phase 
mixing at the origin of the decrease of v is very efficient. In the situation considered by 
O'Neil, electrons with the same initial velocity are not all trapped by the wave, depending 
on their initial position. Moreover, by the time the wave amplitude has reached its constant 
value, the electrons orbits are essentially unperturbed, so that it takes more time for phase 
mixing to be effective. Hence, v is less efficiently reduced to in the O'Neil situation than 
in ours, and we find v « whenever J lo^di > 6, instead of u^t > 30 as found by O'Neil. 
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B. Three dimensional (3-D) space variation 

We now discuss how, and when, 3-D effects may change the results derived previously, 
in the limit of a nearly unperturbed transverse electron motion. In case of a laser driven 
plasma wave, and when the laser electric field is polarized along the y direction, one easily 
finds from Newton equations, 

v y = v 0y + 0(eA/m), (37) 
v z = v 0z + O[(eA/m) 2 /c], (38) 

where A is the amplitude of the laser vector potential, while v 0y and v 0z are the unperturbed 
transverse velocities. Hence, the transverse motion may be considered as unperturbed pro- 
vided that eA/m <C v t h- This condition is fulfilled, for example, for typical laser and plasma 
conditions met in inertial confinement fusion. 

Let us now consider electrons with the same transverse velocities. Their contribution to 
Im(x), which we denote by Iid{vq V) vq z ), is derived from the formulas of Sections!]!! andHVl 
provided that all quantities such as J^UBdt, or 7, be now calculated in the frame moving at 
velocity v = v^x + vo y y + v Qz z with respect to the laboratory frame since, in this frame, the 
electrons have no transverse motion. In particular, J^usdu now is, f Q Ub[x — J v ( / > (t')dt' , y — 
v 0y (t — u),z — v 0z (t — u),u]du, and clearly assumes lower values than in 1-D. Indeed, the 
electrons interact with the wave during a smaller time since, due to their transverse motion, 
they escape more rapidly from the region where the wave amplitude is significant. We 
therefore expect Im(%) to remain close to its linear value, and v close to Ul, up to longer 
times in 3-D than in 1-D. Now, in order to calculate Im(x), we just need to sum over all 
contributions IiD(v 0y ,v 0z ), that is, 

/+00 
Imivoy, v 0z )f (v 0y , v 0z )dv 0y dv 0z , (39) 
-00 

where fo(vo y , Vo z ) is the unperturbed transverse distribution function. Im(%) assumes values 
significantly different from those derived in 1-D if \n y)Z Vth\ ^ |T + nv t h\, where n y ^ z = 
E~ l d y;Z E p , that is when the field amplitude variations experienced by the electrons is mainly 
due to the y or z dependence of E p . Then, not only would v decrease later as a function of 
time, but also more smoothly because the Heavyside-like function found in Sections [III and 
ITVl is now convoluated with / . Hence, v becomes a complicated operator of the transverse 
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gradients of the wave amplitude, and may only be seen again as a damping rate if these 
gradients may be viewed as given parameters. For example, in case of a laser-driven plasma 
wave, the transverse dependence of E p is directly related to that of the laser intensity, due 
to its focusing inside of the plasma, and may therefore be considered as given. 



V. NONLINEAR FREQUENCY SHIFT OF A DRIVEN PLASMA WAVE 



In this Section, we briefly recall results discribed in Ref. |13j regarding the nonlinear 
frequency of a driven plasma wave. Plugging the definition ([7]) of % into Gauss law, one 
finds the following dispersion relation, 

1 + a d Re( X ) = 0, (40) 

where 

l + 2(E d /E p )sm(5<p) + (E d /E p ) 2 
ad l + (E d /E p )sm(8 V ) • 1 } 

When the plasma wave is not driven, and E d = 0, a d = 1 and one recovers the usual 
dispersion relation 1 + Re(x) = 0. The linear value, «ii n , of a d is chosen so as to correspond 
to the linearly most unstable wave against SRS, and its value results from the optimizing of 
two opposite trends. On one hand, it seems clear that it is easier to drive an electrostatic 
wave if this wave is a natural plasma mode. Hence, an n should be close to unity. On the 
other hand, a wave grows more effectively if its Landau damping rate is small, that is if its 
phase velocity is large compared to the thermal one. Since, for a given wave number, k, the 
frequency u derived from Eq. ( 1401) increases with a d , we conclude that an n > 1. Moreover, 
because the Landau damping rate increases with k\n, so does an n . Now, from Eq. (CQ) it is 
clear that, due to the decrease of Im(x) shown in the previous Sections, E d /E p decreases as 
the plasma wave grows, which entails a rapid drop towards unity of a d and hence a rapid 
initial decrease of u. As a consequence, the frequency shift, 5u = u — c^v, where u\ in is 
the EPW linear frequency, is larger in magnitude than could be found by assuming that the 
wave was freely propagating i.e., by solving Eq. ( T4U1) with a d = 1. This is illustrated in 
Fig. [9] which clearly shows that the initial drop in 5u is missed if one assumes a d = 1 when 
solving Eq. (1401) . How to accurately calculate the nonlinear values of a d is explained in Ref. 
[13I and, accounting for the decrease of a d allowed us to derive values of 5uj in very good 
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FIG. 9: The nonlinear frequency shift of the plasma wave, Su, as calculated numerically from 
Vlasov simulations of SRS (green dots), theoretically by solving Eq. (jlQj) (black solid line), and by 
solving Eq. (|40p with ad = 1 (blue dashed line), when kXp ~ 0.52. 

agreement with those derived from Vlasov simulations of SRS, as shown in Fig. [9] when 
kX D « 0.52. 

After the initial drop in u due to that of ct^, the plasma wave frequency keeps on decreas- 
ing due to the nonlinear change in Re(x), which is calculated by making use of the adiabatic 
approximation. Then, the value we find for Re(x) in the limit of a vanishing wave amplitude 
is the same as that published, for example, in Refs. 17, 18, Q, [2^]. However, unlike in these 



papers, we do find solutions to the dispersion relation when kXo > 0.53, and for an infinitely 
small wave amplitude, because we solve Eq. (14"01 with ay ^ 1. Physically this means that, 
by sending a laser into a plasma it is always possible to drive an electrostatic wave, even with 
kXp > 0.53 and slowly enough for an adiabatic estimate of Re(%) to be valid, as shown in 
Ref. [13j. In order to calculate the nonlinear values of Re(x), by making use of the adiabatic 
approximation, we account for the nonlinear change of the phase velocity, which allows us 
to find solutions to the dispersion relation Eq. (1401) up to much larger values than if we had 



assumed that the wave frame was inertial, as was done in Refs. 



18 



20]. 
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VI. APPLICATION TO STIMULATED RAMAN SCATTERING 



In this Section, we briefly discuss how our theoretical model applies to the studying of 
stimulated Raman scattering in the nonlinear regime, and we actually focus on the threshold 



of the so-called "kinetic inflation". This term was used in Ref. [lOj to design the regime 
where SRS reflectivity was experimentally found to be much larger than could be inferred 
from linear theory, a result which was attributed to the nonlinear reduction of the Landau 
damping rate. 

In its simplest version, SRS is a three wave process, an incident laser generating an 
electron plasma wave and a scattered electromagnetic wave. We assume that each of the 
electric field of these waves writes in terms of a slowly varying amplitude and an eikonal 
i.e., that the total electric field is, 

E tot = E p sm(ip p )x + y [E t sin(y?2) + E s cos(y2 s )] , (42) 

where E p , E\ and E s are, respectively, the plasma, laser and scattered wave amplitude. We 
moreover require \E~^ s d t E p ^ s \ < \d t ip p ^ s \ and \E~l s d x E pAs \ < \d x ip p j }S \. Then, in order to 
address the issue of SRS, one actually needs to solve three coupled envelope equations, one 
for each wave. It is actually more convenient to write these equations on complex quantities, 
which lets us define, 

E p = 2E 0p , (43) 
E x = 2E 0l e i(k ' mx -^ nt) e-^', (44) 
E s = 2E 0s e i{kl ™ x -^ nt) e- iVa e i ti 5 " {x ' u)du ) (45) 

where k\ m and k l ™ are the linear values of the laser and scattered wave numbers, k^ s = d x (fi jS , 
uj\ m and u l s in are the linear values of the laser and scattered frequencies, uo^ s = —d t ipi^, and 
8uj is the nonlinear frequency shift of the plasma wave, defined in Section |V] Using Maxwell 
equations, and writing Gauss law as described in the previous Sections, we find the following 
equations, valid for a uniform plasma and in 1-D, 

dE Qp i „. dE Qp i „p Re(T p E iE£ s ) 

~dT + V ^ + VE ^ ~ dujXlnv ' (46) 

^ 1 + v gs ^^ + i[5uj-Vg S 5k]EQ S = T s E 01 Eq p , (47) 

it +v "it = - v > E °- E °" < 48 ' 
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FIG. 10: Reflectivity, R, as a function of time when the laser intensity is // = 1.5 x 10 14 W/cm 2 , 
panel (a), as calculated using the Vlasov code ELVIS and, panel (b), using the envelope code 
BRAMA, and when the laser intensity is I\ = 2 x 10 14 W/cm 2 , panel (c), as given by a Valsov 
simulation and, panel (d), as given by our envelope code. 



where, in Eq. (1471) . 5k is the nonlinear wave number shift of the plasma wave, related to Sco 
by the equation, d t Sk = —d x 5ui, v g i and v gs are the usual group velocities of electromagnetic 
waves, v gs = kic/ui, v gs = k s c/u s , as for T p , T/ and T s , these are plain constants, T p = 
ek/mujiuj s , T s = ek/2muJi, and Ti = ek/2muj s , where k = d x (p p is the plasma wave number. 
The envelope equations (I46H48P are solved using the code BRAMA, which will be detailed 
in a forthcoming paper, and the results are compared to those of the Vlasov code ELVIS, 



Ref. 191. In our simulations, either with the Vlasov or the envelope code, SRS results from 
the optical mixing of a laser, and a counterpropagating seed, as explained in Section HVl 
The ratio between the seed intensity I S {L) at the right end of the simulations box, and 
the laser intensity at the left end of the box, Il(0), is chosen to be 10~ 5 . Figure [TUI plots 
the reflectivity R = I s (0)/ Il{0) as a function of time, calculated for a 1-D uniform plasma 
with electron temperature, T e = 2keV, electron density n = 0.1n c , and whose length is 
lOOyum. The laser wavelength is 0.35 /im while the seed wavelength is 0.55 /im. When the 
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laser intensity is J/ = 1.5 x 10 14 W/cm 2 , a linear theory would predict the reflectivity to be 
Ru n pa 2 x 10~ 5 , and both the Vlasov and envelope codes find R of the order of 10~ 5 . By 
contrast, when l\ = 2 x 10 14 W/ cm 2 , while the linear value of the reflectivity is Rn n pa 3 x 10~ 5 , 
the reflectivity calculated either with the Vlasov or the envelope code is of the order of 10%, 
as can be seen in Fig. [TUJ This Figure also shows some discrepancies in the actual values 
of the reflectivity predicted by the two different codes, whose origin will be discussed in a 
future paper and is way beyond the scope of this article. However, as regards the threshold 
for inflation, both codes agree that the threshold intensity lies between 1.5 x 10 14 W/cm 2 
and 2 x 10 14 W/cm 2 , while the envelope code is about 5000 faster in providing this result. 
Hence, using the theoretical model described in the previous Sections, we built a powerful 
tool to predict when stimulated Raman scattering is negligible, which an important issue 
for inertial confinement fusion (see for example Ref. [3]). 



VII. CONCLUSION 



In this paper, we investigated how efficiently an electron plasma wave (EPW) could be 
externally driven. This led us define the nonlinear group velocity, v g , and Landau damping 
rate, z/, of a driven plasma wave, which are terms appearing naturally in the envelope 
equation for the wave amplitude. We provided a practical analytic formula for z/, and found 
the unexpected result that v assumed nearly constant values before abruptly dropping to 
zero, and that this drop in v occurred simultaneously with a rapid increase of v g towards 
the wave phase velocity, and a decrease of the coupling constant between the plasma wave 
and the driving field. We moreover unambiguously showed, without resorting to complex 
contour deformation, that a plasma wave, first driven by laser at a small enough amplitude 
and then freely propagating, would damp at the rate predicted by Landau. This then 
imposes restrictions for non-Landau damping, as predicted by Belmont et al. in Ref. {fl, to 
indeed occur in actual experiments. All these results stem from our theoretical derivation of 
Im(x), which directly follows from the investigation of the nonlinear electron motion. The 
expression found for Im(x) actually results from the matching of two very different estimates, 
a perturbative one for small amplitudes, and one relying on the adiabatic approximation 
and valid whenever v pa 0. This yields values for Im(x) in excellent agreement with those 
either inferred from test particle simulations or from Vlasov simulations of stimulated Raman 
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scattering (SRS). 

We moreover discussed in this article the nonlinear frequency shift, 8u, of a driven plasma 
wave and found that \Su\ was much larger than could be derived by assuming that the wave 
was freely propagating. We moreover showed that no physical effect could be attributed 
to the increase of k\p above 0.53, unlike what could be inferred from Ref. This 
emphasizes the importance of specifying the way a plasma wave has actually been generated 
in order to discuss its nonlinear properties. 

Our results regarding both, the EPW envelope equation and its nonlinear frequency 
shift, allow us to study SRS in the nonlinear regime. In particular, we investigated the 
threshold of the so-called kinetic inflation, a regime where the SRS reflectivity is much 
larger than predicted by linear theory. This threshold is a very important parameter for 
inertial confinement fusion because, below it, one is assured that SRS reflectivity would be 
very low and therefore that SRS would not affect the fusion efficiency. Using our model 
when the plasma is homogeneous, and in a 1-D geometry, we found values for the inflation 
threshold in very good agreement with those derived from Valsov simulations, but within 
a much smaller computing time. This shows the potentiality of our model to address more 
complicated physics situations. 

In conclusion, we derived very precisely the nonlinear properties of a driven electron 
plasma wave, which allowed us to discuss the generality of previous results on this topic, 
which is a long standing, and basic issue in plasma physics. We moreover applied our results 
to the studying of stimulated Raman scattering, and to the threshold for kinetic inflation, 
which is an important issue for inertial confinement fusion. 

APPENDIX A: HAMILTON! AN PERTURB ATIVE ANALYSIS 

In this Appendix, we use a first order Hamiltonian perturbative analysis to approximate 
the motion of an electron acted upon by a longitudinal wave whose electric field is E = 
£7o(^) el1 ^ +c.c, and whose frequency, to, and wave number, k, are defined by k = d x p, to = 
—d t (p. In the dimensionless variables, r = t/kv t h, <p{t) = (p[x(r), r] and v = v^dx/dt, where 
Vth = y/T e /m is the thermal velocity, the electron dynamics derives from the Hamiltonian, 




(Al) 
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where $ = eE§jkT ei and = uo/kvth- The perturbative calculation consists in defining a 
canonical change of variables (ip, v) — > (ip', v') such that v' is a constant of motion, at least 
at first order in the wave amplitude. The change of coordinates is defined using a generative 
function, F(ip,v'), and is 

ip' = ip + d v ,F, (A2) 
v = v' + d v F. (A3) 

Then, ip « ip + (vq — v^r + Sip, where ip and v are constant, and 

Sip = -d' v F. (A4) 

In the new variables, the new Hamiltonian is, 

ot 2 ot 

The generative function, F, is then chosen so as to cancel the term i$>e ltp + c.c., so that, at 
first order in $, it needs to solve, 

dF OF 

(v' - v*) — + — = + c.c. (A6) 

dip ot 

We now assume that, at r = 0, the wave amplitude is infinitesimal, so that Sip = F = 0. 
Then, the solution of Eq. ( 1A6I) is, 



F = -ie i(p / $(M)e™ (u - r) .c.c, (A7) 
Jo 

where we have denoted w = v' — v$. Then, 

5ip = -ie^d w (J $(M)e™ (u - r) .c.c.^ 

« -ie i{lp0+WT) d w Qf <$>(u)e lw(u - T) .c.cA (A8) 

APPENDIX B: APPROXIMATE EXPRESSION FOR Im( Xa )- 

In this Appendix, we give an approximate expression of 



Xa 



(fcA D ) 2 $(r) 



-r) / iwe iu,(f - T) «M, (B.l) 

«/|«u|>Vj 
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in the limit Vj C r^ 1 , where r$ is the typical timescale of variation of $. From the results 
of Section [TTTj. it is clear that Im(xa) = — (k\o) ~ 2 fo{ v ^)[^ + 8Xa], with, 

5 Xa = Hr)- 1 f (u- r)®{u)d u G{u - r)du (B.2) 
Jo 

where 

rVi 

G(u -t)= e iw(u - T) dw. (B.3) 

J-Vi 

Clearly, the timescale of variation of G is Vf 1 , while d u G\ u=T = 0, and d 2 2 G\ u=T = 2VJ 3 /3. 
Then, integrating (1B.2|) three times by parts yields 

^3 FT FU {■£ 



4V, 6 r f u n 

+<S>(t)-> f [(« - T)aJ,G + ( f" t f (B.4) 

Jo Vjo Jo Jo I 



Clearly, the last term in the right-hand side of Eq. flB.41) is of the order (VJr$) times the 
first one, and is therefore negligible in the limit V\ r^ 1 . 
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